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Abstract 



Evaluating the neutronic state of the whole nuclear core is a very im- 
portant topic that have strong implication for nuclear core management 
and for security monitoring. The core state is evaluated using measure- 
ments. Usually, part of the measurements are used, and only one kind 
of instruments are taken into account. However, the core state evalua- 
tion should be more accurate when more measurements are collected in 
the core. But using information from heterogeneous sources is at glance 
a difficult task. This difficulty can be overcome by Data Assimilation 
techniques. Such a method allows to combine in a coherent framework 
the information coming from model and the one coming from various 
type of observations. Beyond the inner advantage to use heterogeneous 
instruments, this leads to obtain a significant increasing of the quality 
of neutronic global state reconstruction with respect to individual use of 
measures. In order to present this approach, we will introduce here the 
basic principles of data assimilation focusing on BLUE (Best Unbiased 
Linear Estimation). Then we will present the configuration of the method 
within the nuclear core problematic. Finally, we will present the results 
obtained on nuclear measurement coming from various instruments. 

Keywords: Data Assimilation, neutronic, activities reconstruction, 
nuclear in-core measurements 

1 Introduction 

The knowledge of the neutronic state in the core is a fundamental point for the 
design, the safety and the production processes of nuclear reactors. Due to the 
crucial role of this information, considerable work has been conducted for long 
time to accurately estimate the neutronic spatial fields. Spatial distribution of 
power or activity in the whole core, or hottest point of the core, can be derived 
from such spatial fields. These information allows mainly to check that the 
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nuclear reactor is working as expected in a very detailed manner, and that it 
will remain in the operating limits when in production. 

Two types of information can be used for the neutronic state evaluation. 

First, the physical core specifications, including the nuclear fuel description, 
make possible to build a numerical simulation of the system. Taking into account 
neutronic, thermic and hydraulic spatial properties of the nuclear core, such 
numeric models calculate in particular the useful reaction rates for the physical 
analysis of the core state. 

Secondly, various measurements can be obtained from in-core or out-of-core 
detectors. Some detectors can measure neutron density, either locally or in 
spatially integrated areas, other can measure temperature of the in-core water 
at some points, etc. A lot of reliable measures comes from periodical flux maps 
measured in each core reactor, at a periodicity of about one month. Then, 
all these measurements do not have the same type of physical relation with 
the neutronic activity, and also not the same accuracy. So it is uneasy to 
take into account simultaneously all these heterogeneous measurements for the 
experimental evaluation of the neutronic state in the core. 

A lot of these measures are local, in determined fuel assemblies, and do 
not give informations in un-instrumented areas of the core. The activity dis- 
tribution over the whole core is traditionally obtained through an interpolation 
procedure, using the calculated field as first guess (a proxy) of the "real" activ- 
ity field corresponding to the measurements. In other words, the activity value 
in un-instrumented areas is calculated as the weighted average of the activity 
measures, using the calculated activity field to interpolate. The power is then 
obtained from activity through an observation operator, which depends only on 
core nominal physical specifications for the periodical flux map measurements. 
This interpolation procedure gives already good results, but some drawbacks 
remains in using only activity measures in a deterministic interpolation proce- 
dure. 

Both physical core specifications and real measurements are subject to some 
uncertainties. Moreover, numerical assumptions, required to use the models, 
add some inaccuracy. All these uncertainties are not used explicitly in the in- 
terpolation procedure, but often used to qualify the a posteriori activity field 
build by the procedure. Moreover, the interpolation cannot take into account, 
for example, heterogeneous instruments, or observed discrepancy of some in- 
struments. 

Attempts have been made to overcome these limitations, mainly in two direc- 
tions. First, studies attempt to combine activity measurements and calculations 
through least-squares derived methods (for example in pQ), leading to the most 
probable activity (or power) field on the whole core. These methods allow to 
take into account heterogeneous measures, but are difficult to develop because 
of their extreme sensitivity to the weighting factor in the combination of mea- 
sures and calculations. Secondly, explicit control of the error, in order to reduce 
its importance, have been tried through the development of adaptive methods 
to adjust coefficients in the calculation or the interpolation procedure. 

Some of these difficulties can be solved by using data assimilation. This 
mathematical and numerical framework allows combining, in an optimal and 
consistent way, values obtained both from experimental measures and from a 
priori models, including information about their uncertainties. Commonly used 
in earth sciences as meteorology or oceanography, data assimilation has strong 
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links with inverse problems or Bayesian estimation [H [3] . It is specifically tai- 
lored to solve such estimation problems through efficient yet powerful procedures 
such as Kalman filtering or variational assimilation |4 j. Already introduced in 
nuclear field 0|6l[7], it can be used both for field reconstruction or for parameter 
estimation in a unified formalism. 

Data assimilation can treat information coming from any type of measure in- 
struments, taking into account the way the measure is related with the objective 
field to be reconstructed, such as neutronic activity here. Data assimilation can 
further adapt itself to instrument configuration changes, and for example the 
removal or the failure of an instrument. Moreover, the method takes natively 
into account informations on instrumental or model uncertrainties, introduc- 
ing them a priori through the data assimilation procedure, and obtaining a 
posteriori the reduced uncertrainties on the reconstruction solution. 

In this paper, we introduce the data assimilation method and how it address 
physical field reconstruction. Then we will make a detailed description of the 
various components that are used in data assimilation, and of the various type of 
instruments we can use to make in-core neutronic activity measurement. Then 
we present results with various instrument situations in nuclear core, obtained 
on a set of true nuclear cores. 

2 Data assimilation 

We briefly introduce the useful data assimilation key points, to understand 
their use as applied here. But data assimilation is a wider domain, and these 
techniques are for example the keys of nowadays meteorological operational 
forecast. This is through advanced data assimilation methods that long-term 
forecasting of the weather has been drastically improved in the last 30 years. 
Forecasting is based on all the available data, such as ground and satellite 
measurements, as well as sophisticated numerical models. Some interesting 
information on these approaches can be found in the following basic references: 

The ultimate goal of data assimilation methods is to be able to provide a best 
estimate of the inaccessible "true" value of the system state, denoted x*, with 
the t index standing for "true" . The basic idea of data assimilation is to put 
together information coming from an a priori state of the system (usually called 
the "background" and denoted x b ), and information coming from measurements 
(denoted as y). The result of data assimilation is called the analysed state x a 
(or the "analysis"), and it is an estimation of the true state x* we want to find. 
Details on the method can be found in [3] or [3]. 

Mathematical relations between all these states need to be defined. As the 
mathematical spaces of the background and of the observations are not necessary 
the same, a bridge between them has to be build. This bridge is called the 
observation operator H, with its linearisation H, that transform values from 
the space of the background state to the space of observations. The reciprocal 
operator is known as the adjoint of H. In the linear case, the adjoint operator 
is the transpose H T of H. 

Two additional pieces of information are needed. The first one is the relation- 
ships between observation errors in all the measured points. They are described 
by the covariance matrix R of observation errors e Q , defined by e G = y — H (x*). 
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It is assumed that the errors are unbiased, so that E[e Q ] = 0, where E is the 
mathematical expectation. R can be obtained from the known errors on the 
unbiased measurements. The second one is similar with the relationships be- 
tween background errors. They are described by the covariance matrix B of 
background errors e&, defined by = x fc — x*. This represents the a priori 
error, assuming it to be unbiased. There are many ways to obtain this a priori 
and background error matrices. However, they are commonly build from the 
output of a model with an evaluation of its accuracy, and/or the result of expert 
knowledge. 

It can be proved, within this framework, that the analysis x a is the Best 
Linear Unbiased Estimator (BLUE), and is given by the following formula: 

x a =x fc + K(y-Hx h ) (1) 

where K is the gain matrix |J: 

K = BH T (HBH T + R) 1 (2) 

Moreover, we can get the analysis error covariance matrix A, characterising the 
analysis errors e Q = x a — x*. This matrix can be expressed from K as: 

A = (I - KH)B (3) 

where I is the identity matrix. 

The detailed demonstrations of those formulas can be found in particular in 
the reference j4]. We note that, in the case of Gaussian distribution probabilities 
for the variables, solving equation [l] is equivalent to minimising the following 
function J(x), x a being the optimal solution: 

J(x) = (x-x b ) T B- 1 (x-x fc ) 

T (4) 
+ (y-Hx^R-^y-Hx) 

We can make some enlightening comments concerning this equation [4j and 
more generally on the data assimilation methodology. If we do extreme assump- 
tions on model and measurements, we notice that these cases are covered by 
minimising J. First, assuming that the model is completely wrong, then the 
covariance matrix B is oo (or equivalently B _1 is 0). The minimum of J is then 
given by x a = H x y (denoting by HU 1 the inverse of H in the least square 
sense). It corresponds directly to information given only by measurements in 
order reconstruct the physical field. Secondly, on the opposite side, the assump- 
tion that measurements are useless implies that R is oo. The minimum of J 
is then evident: x a = x b and the best estimate of the physical field is then de 
calculated one. Thus, such an approach covers whole range of assumptions we 
can have with respect to models and measurements. 

3 Data assimilation method parameters 

The framework of the study is the standard configuration of a 1300 electri- 
cal MW Pressurized Water Reactor (PWR1300) nuclear core. Our goal is to 
reconstruct the neutronic fields, such as the activity, in the active part of the 
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nuclear core. For that purpose, wc use data assimilation. To implement such 
methodology, we need both simulation codes and measures. For the simulation 
code, we use standard EDF calculation code for nuclear core simulation, in a typ- 
ical configuration. The results are build on a set of 20 experimental neutronic 
flux maps measured on various PWR1300 nuclear cores. Such measurements 
are done periodically (about each month) on each nuclear core. These different 
measurement situations are chosen for their representativeness, in order that 
statistical results covers a wide range of situations and can have some sort of 
predictability property. 

3.1 The backgroud and the measurements 

A standard PWR1300 nuclear core has 193 fuel assemblies within. For the cal- 
culation, those assemblies are each considered as homogeneous, and are divided 
in 38 vertical levels. Thus, the state field x can be represented as a vector of 
size 193 x 38 = 7334. 

The measurements come from instruments that can be located on horizontal 
2D maps of the core. There are three types of instrument that are usually used 
to monitor the nuclear power core: 

• Mobile Fission Chambers (MFC), which are inside the active nuclear core, 
for which the locations are presented on the figure [T] 

• Thermocouples (TC), which are above the active nuclear core, for which 
the locations are presented on figure [2] 

• fixed ex-core detector locations. 

The data coming from the ex-core detectors are continuous in time and are 
very efficient for security purpose, which is their main goal. Their purpose is 
to continuously monitor the core, but not to measure accurately the neutronic 
activity at each fine flux map. So, their measures are too crude for being 
interesting on a fine reconstruction of the inner core activity map. Thus, we do 
not take into account information coming from those data. 

All these types of instrumentation (MFC, TC, ex-core) can be found on any 
power plants. For the purpose of this study, we add artificially an extra type 
of detector, described as idealized Low Granularity MFC (named here LMFC). 
The measurement attributed to the LMFC are build artificially from the in- 
formation given by the MFC. Thus they are replacing the MFC on the given 
LMFC locations. The evaluation of LMFC response is calculated from the MFC 
measured neutron flux, assuming a different physic process, and a lower granu- 
larity. The lower granularity assumption done on the LMFC induces a partial 
integration of the results of the MFC over a given area. Of course, the phys- 
ical process involved to make a measurement being different, the resolution of 
LMFC will be different from the one from MFC. We take 16 of those instru- 
ments. They are located in various area of the core, replacing MFC, to try to 
make a representative array of measurement as in is shown on figure [3] 

The main characteristics of the instruments, as there number in the core, the 
number of considered vertical levels, and the size of the part of the observation 
vector y associated with the particular instrument type, are reported in table 
[l] The size of the final observation vector is given by summing the size of all 
the individual y vector of the instruments used. 
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Figure 1: The Mobile Fission Chambers (MFC) instruments within the nuclear 
core are localised in assemblies in black, in a horizontal slice of the core. The 
assemblies without instrument are marked in white, and the reflector is in grey. 



Instrument 


Locations 


Vertical 


Size of the y 


type 


number 


levels 


vector part 


MFC 


42 


38 


1596 


LMFC 


16 


8 


128 


TC 


50 


1 


50 



Table 1: Main characteristics of the instruments used for data assimilation. 
These characteristics remains the same, either in mono-instrumented cases or 
in multi-instrumented ones. 



3.2 The observation operator H 

The H observation operator is the composition of a selection and of a normali- 
sation procedure, and can be build independently for each instrument. So there 
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Figure 2: The Thermocouples (TC) instruments within the nuclear core are 
localised above assemblies in black, in a horizontal slice of the core. The as- 
semblies without instrument are marked in white, and the reflector is in grey. 



is one individual H matrix observation operator by instrument type. The com- 
plete H matrix observation operator is the concatenation, as a bloc-diagonal 
matrix, of all the individual matrix for each instrument. 

Each observation operator is basically a selection matrix, that choose in the 
model space a cell that is involved in a measurement in the observation space. 
In addition, a weight, according to the size of the cell, is affected to the selection. 
As experimental data are normalised, this selection matrix is multiplied by a 
normalisation matrix that represents the effect of the cross normalisation of the 
data. This observation matrix is a (7334 x dj) matrix, where di is the size of 
the part of the observation vector y for each instrument involved in assimilation, 
as reported in table [T] 
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Figure 3: The idealized Low Granularity MFC (LMFC) instruments within the 
nuclear core are localised in assemblies in black, in a horizontal slice of the core. 
The assemblies without instrument are marked in white, and the reflector is in 
grey. 



3.3 The background error covariance matrix B 

The B matrix represents the covariance between the spatialised errors for the 
background. The B matrix is estimated as the double-product of a correlation 
matrix C by a diagonal scaling matrix containing standard deviation, to set 
variances. 

The correlation C matrix is built using a positive function that defines the 
correlations between instruments with respect to a pseudo-distance in model 
space. Positive functions allow (via Bochncr theorem) to build symmetric de- 
fined positive matrix when they are used as matrix generator (for theoretical 
insight, see reference documents [TU] and dJ). Second Order Auto-Regressive 
(SOAR) function is used here. In such a function, the amount of correlation 
depends from the euclidean distance between spatial points in the core. The ra- 
dial and vertical correlation lengths (denoted L r and L z respectively, associated 
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to the radial r coordinate and the vertical z coordinate) have different values, 
which means we are dealing with a global pseudo euclidean distance. The used 
function can be expressed as follow: 

The matrix C obtained from the above Equation[5]is a correlation one. It can be 
multiplied (on left and right) by a suitable diagonal standard deviation matrix, 
to get covariance matrix. If the error variance is spatially constant, there is only 
one coefficient to multiply C. This coefficient is obtained here by a statistical 
study of difference between the model and the measurements in real case. In 
real cases, this value is set around a few percent. 

The size of the background error covariance matrix B is related to the size 
of model space, so it is (7334 x 7334) here. 

3.4 The observation error covariance matrix R 

The observation error covariance matrix R is approximated by a simple diagonal 
matrix. It means we assume that no significant correlation exists between the 
measurement errors of all the instruments. A usual modelling consists in taking 
the diagonal values as a percentage of the observation values. This can be 
expressed as: 

R ij = {°Vjf > Vj. (6) 

The a parameter is fixed according to the accuracy of the measurement and 
the representative error associated to the instrument. It is the same for all the 
diagonal coefficients related with one instrument. Its values is only depending 
of the type of instrument we are dealing with. The a value can be determined 
by both statistical method and expert opinion about the measurement quality. 
In the present paper, we will use arbitrary value for the a. 

The size of the R matrix is related to the size of the observation space, 
so it is (J^i di x Y^,i di) where di is the size of the observation vector of each 
instrument i involved in assimilation, as reported in the table [T] 

4 Results on data assimilation using only one 
type of instrument 

The first results are showing the quality of the reconstruction as a function of 
the various types of instruments that are taken into account for reconstruct the 
activity of the core. 

The experimental data are a set of measurement on the 38 levels of the all 
the instrument locations inside of the core. Thus, to evaluate the quality of the 
reconstruction of the physical fields with one type of instrument, we look for the 
misfit (y — Hx a ) at measurement locations (by other instruments) that are not 
involved in the assimilation process. The number of locations, where there is a 
measurement and that is not involved in data assimilation procedure, so where 
the misfit (y — Hx") is calculated, are synthesised in the table [2j as well as the 
accuracy associated to each instrument through the a parameter of equation [6] 
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Instrument 
type 


Number of misfit 
calculations locations 


a value 


MFC 


16 


1% 


LMFC 


42 


2% 


TC 


58 


3% 



Table 2: Number of misfit calculation locations used for each instrument type 
considered individually for data assimilation, and arbitrary accuracy assumed 
in the present studies for each type of instrument. 

Thermocouples being a fully integral measurement outside of the active core, 
we can do the misfit (y — Hx") evaluation of the reconstruction in all the 
locations of MFC/LMFC. 

For each data assimilation procedure associated with an instrument, we cal- 
culate the Root Mean Square (RMS, which is the norm) of the misfit (y — Hx") 
on all the misfit calculation locations. To synthesize the value for one set of 
measurement, we take the mean value of the misfit (y — Hx") on each of the 
38 levels of the core, which leads to a horizontally averaged value of the misfit 
(y — Hx"). In sake of more general behaviour, we take 20 set of flux map mea- 
surements, with various settings and ageing of PWR1300 nuclear cores. Then 
we proceed to the calculation of the mean value on all those set of measurement. 

Globally all the results present strong misfit (y — Hx") on the upper and 
lower levels, which is a known effect mainly due to the axial reflector modelling. 
In addition, the centre part of the reactive core in the nuclear plant is also 
the region where the neutron flux is the most intense, so the hot spot of the 
activity field in the core is expected to be in this region. Thus, the next plots 
are restricted to the centre of the core. 

Figure [3] shows the axial misfit measured by the standard RMS of the dif- 
ference between analysis and measurements, in arbitrary units, on all the data 
assimilation unused locations for the various types of instrument studied. The 
oscillating behaviours, that are barely noticeable on all the curves, come from 
the different material that are within a core level. Some levels are containing 
mechanical structure of the core, thus these are more neutron absorptive. 

We noticed, as expected, that the reconstruction coming from the thermo- 
couples (TC) is the closest to the background, due to their integral measurement 
property and their lower accuracy. Moreover, an improvement of the accuracy 
does not improve dramatically the quality of the core state evaluation, mainly 
due to the integral measurement property. 

The Low Granularity MFC (LMFC) are showing a good reconstruction of 
the physical field in the whole core, despite the not so good accuracy and the 
limited number of measurement. The increase of the misfit from (y — Hx 6 ) to 
(y — Hx") for the lower part of the core is easily explainable by the chosen 
locations of LMFC that are in core. This part, near the border of the core, does 
not get enough measurements to be very accurately reconstructed. 

The reconstruction using only MFC is, also as expected from accuracy values, 
the best one. The data assimilation procedure leads to half the misfit (y — Hx 6 ) 
observed when only using the model. 

From results coming from MFC and LMFC, we notice that, within the hy- 
pothesis and the chosen modelling of the integration, TC measurement are per- 
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Figure 4: Vertical misfit for various kind of intruments, measured by the RMS of 
the differences between the analysis and the measurements at unused locations 
(see text for details). 



mitting only a crude evaluation of the core state. 

5 Results on data assimilation with heteroge- 
neous instruments 

This section describes results using different instruments together in the data 
assimilation procedure. 

In this case, we cannot take as reference the measures points taken apart, 
as when we are studding individual instrument as presented in section [4] Thus, 
we choose to make an analysis with data assimilation using all the available 
measures in the core, on the 58 locations, and using a 1% accuracy. Then 
we evaluate the misfit (y — Hx a ) with respect to this reference calculation in 
all the assemblies where no MFC are present, which means 135 locations (193 
total assembly minus 58 instrumented locations) . The calculated values at those 
instrumented locations allow to benchmark the quality of the reconstruction. On 
this misfit (y — Hx"), we calculate the RMS per horizontal level as previously, 
and take the average over some 20 selected flux map measurements. The results 
are presented on figure [5j 

The figure [5] presents RMS per horizontal level of the misfit calculated for 
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Figure 5: Axial misfit for various kind of instrumental configurations, measured 
by the RMS of the difference between the analysis and the reference calculation. 



various instrument taken alone or in conjunction. We notice that the results 
have the same behaviour as the one we got in hgure |4j This means that the 
reference we chose to evaluate the quality of data assimilation using several 
instrument types can be considered as reliable. 

Looking at the successive addition of instrument, we notice that addition of 
the LMFC have a important effect, the MFC+LMFC configuration presenting 
an improvement with respect to the configuration with only MFC. However, 
adding thermocouple to this configuration is not really helpful, and the im- 
provement is not really noticeable on this figure [5] 

These results highlight a number of important points on data assimilation 
methodology. On one hand, when only very few measurements are available, 
they are very helpful and allow a fairly good improvement of reconstructed 
physical field. On the other hand, when a lot of measurements are available, 
adding a few more, or with a lower accuracy as thermocouples, do not change 
dramatically the result of field reconstruction by data assimilation. On overall, 
this shows that data assimilation technique is doing the best use of experimental 
information provided to the procedure. Those results are comforting the ones 
found on the robustness of the evaluation of the nuclear core by data assimi- 
lation, when only MFC are used as presented in [7j. Moreover, as expected in 
data assimilation technique, the use of heterogeneous instrument is transparent 
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within the method. 



6 Conclusion 

The use of data assimilation has already been proved to be efficient to recon- 
struct fields in several domains, and recently in neutronic activity field recon- 
struction for nuclear core. The present paper demonstrates that, within the data 
assimilation framework, information coming from heterogeneous sources can be 
used in a transparent way without making any adjustment to the method. 

Looking at the various types of instruments, we have (MFC, LMFC and 
TC) we notice that the influence they have on reconstruction depends on three 
parameters: 

• the granularity of each type of instrument, that is the density of instru- 
ments, their integral measurement property and their repartition all over 
the core, 

• the accuracy of each instrument, possibly with respect to the accuracy of 
the others, 

• and the global instrumentation configuration, that is the complex repar- 
tition of all instruments. 

Those conclusion arise from comparison of the various instruments, using them 
individually or together, in various instrumental configurations. 

Data assimilation gives a very efficient and adaptable framework in order to 
take the best from both experimental data and model. Moreover, this can be 
done without heterogeneities constrains between instruments. This technique 
is used here in a very elementary way, but it opens the door to many develop- 
ments, for example in systematic data analysis, models comparison, in dynamic 
modelling, etc. for nuclear reactors. 
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